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Q\ ■ Abstract 

Qh! The direct production rate of J/ip via the color-octet mechanism is calculated at 

the Z resonance. The color-octet production process Z — > J/ipqq is shown to have a 



substantial branching ratio as well as a distinctive energy spectrum, which can be used 



^ ■ as a powerful tool to distinguish from the color-singlet direct production of the J/ip. 

.«?; 

Despite charmonium and bottomonium have been discovered some times ago, heavy 
quarkonium physics has constantly drawn a lot of attentions from theorists. This is mainly 
due to the property of asymptotic freedom of QCD which allows theorists to calculate pertur- 
batively the short- distance part in the production and in the decay rates of quarkonia, while 
the incalculable long-distance part can be factored out and absorbed into the wave-function 
for the ground state or its derivatives for higher orbital excitations. Confinement implies that 
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quarkonium must be seen as a color-singlet object. The most natural assumption to make in 
the calculation of the short- distance part of quarkonium production or decay rate is that the 
heavy quark- ant iquark pair inside the bound state is produced or annihilated with the same 
orbital angular momentum, spin, and color of the asymptotic physical quarkonium states. 
This simple picture dubbed as the color-singlet model in the literature is derivable from 
QCD and has greater predictive power than other alternatives |2| like the color-evaporation 
model or the local duality approach. 

During the past several years, a wealth of quarkonium data collected by the CDF detector 
at the Tevatron has posted great challenge to the naive color-singlet model. With the advance 
of the technology of silicon vertex detector, CDF in the 1992-1993 run was able to separate 
their quarkonium data into two vast classes ||: the first class can all be traced back to 
a secondary displaced vertex which is a few microns away from the interaction point and 
presumably they are all coming indirectly from B meson decays; the second class of the 
data was called prompt because it comes from direct QCD production. Surprisingly, the 
rates for the prompt ip, ip', and Xcj at the large transverse momentum region were observed 
to be orders of magnitude above the lowest order perturbative calculation based on the 
color-singlet model ||. 

To resolve these discrepancies one needs to seek for new production mechanisms as well as 
new physical insights that go beyond the color-singlet model. The factorization model based 
on non-relativistic QCD (NRQCD) recently advocated by Bodwin, Braaten, and Lepage M 
provides a new framework to calculate the inclusive production and decay rates of quarko- 
nium to any order in strong coupling constant a s , as well as to any order in v 2 where v is the 
typical relative velocity of the heavy quarks inside the bound state. At the same time, the 
importance of fragmentation contributions to the quarkonium production at the Tevatron 
was also pointed out by Braaten and Yuan ||. Parton fragmentation into quarkonium is 
formally of higher order in strong coupling constant but they can be enhanced over the tra- 
ditional lowest order gluon-fusion mechanism || at sufficiently high transverse momentum. 
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However, cross sections for prompt ip and ip' production were observed to be more than an 
order of magnitude above the predictions based on the color-singlet model even after the 
inclusion of the fragmentation contributions 0, ||, ||| . 

One crucial feature of the factorization model is that quarkonium is not solely regarded 
as simply a quark-antiquark pair but rather a superposition of Fock states. In the light-cone 
gauge, the Fock state expansion for the ip is 

|V>> = \cc( 3 S 1 ,l)) + 0(v) \cc( 3 Pj,^ + 0(v 2 )(^\cc( 3 S 1 ,8 or l)gg) + \cc( 1 S ,8)g)) 
+ 0(v 2 ) \cc( 3 Dj,,8 or l)gg) + 0(v 3 ) , (1) 

where the angular momenta of the cc pair in each Fock state is labeled by 2S+1 Lj with 
a color configuration of either 8 or L A novel idea proposed by Braaten and Fleming 
[TT| to explain the CDF data on prompt ip (?/>') production is due to a color-octet term 
in the gluon fragmentation function into ip (?//). This term corresponds to the case where 
the fragmenting gluon forms a cc pair in the color-octet 3 S\ state at the short- distance of 
order l/m c or smaller, and subsequently evolves non-perturbatively into the physical ip or 
ip'. The non-perturbative effects are represented by the NRQCD matrix element (0$ ( 3 5'i)) 
((0$ ( 3 Si))) H • These color-octet matrix elements are suppressed by order v A relative to 
the corresponding color-singlet matrix elements (Of ( 3 Si)) ((Of ( 3 Si))). On the other hand, 
the short- distance factor computed in perturbative theory for this color-octet process is only 
of order a s (2m c ), which is enhanced by a factor of 1/atg relative to the color-singlet process. 
Therefore, the suppression in the color-octet matrix element can be easily compensated by 
the enhancement in the corresponding short- distance factor when compared with the color- 
singlet process. The CDF data for ip (ip 1 ) can be explained by including this color-octet 
term in the gluon fragmentation function and by adjusting the value for the matrix element 
(Ot( 3 S 1 )) ((Of( 3 S l ))) to fit the data 0, p|. 

While the color-octet mechanism can easily explain the transverse momentum spectrum 
of the quarkonium measured at the Tevatron, the normalization cannot be determined a 
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priori. One has to rely upon the CDF if> and if)' data to extract the unknown color-octet 
matrix elements. In order to establish firmly our belief of the color-octet mechanism, it 
is therefore important to identify these color-octet quarkonium signals in other production 



processes. Recently, Cho and Leibovich [Tj have studied color-octet quarkonium production 
at the Tevatron using exact tree-level calculations rather than the fragmentation approxi- 
mation. This approach allows them to probe the CDF charmonium and bottomonium data 
at almost the whole range of transverse momentum. Braaten and Chen [H| showed that the 
color-octet mechanism can give rise to a distinctive signature at the upper endpoint of the 
energy distribution of if) at CLEO energy. In this Letter, we study the prompt quarkonium 
production by the color-octet mechanism at the Z resonance. This work was motivated by 
the capability of the LEP detectors to separate the prompt if> from those coming from B 



decays |14|]. Besides, the preliminary results of prompt if) production at LEP |TJ| suggested 
almost a factor of ten larger than the prediction of the most dominant color-singlet process 
via charm quark fragmentation [To*, IE). We shall show that the most important color-octet 



process will give rise to a branching ratio substantially larger than the most dominant color- 



singlet process and is consistent with the preliminary LEP results [13]. Furthermore, the 
energy spectrum predicted by the color-octet production process is very soft and should 
be distinguishable from that produced by the most dominant color-singlet process. Hence, 
our consideration of if> production in Z decays is crucial in verifying the importance of the 
color-octet mechanism in quarkonium production. 

According to the factorization formalism Q], the differential rate for the inclusive produc- 
tion of a charmonium state H with momentum P from Z decay can be written in a factored 
form: 

dT(Z -> H{P) + X)=J2 dT(Z -> cc(P, n) + X) (O*) , (2) 

n 

where dt is the partonic decay rate for producing a cc pair with total momentum P, vanishing 
relative momentum, and in an angular momentum and color state labeled collectively by n. 
The long-distance factor (O^ ) is the NRQCD matrix element describing the probability for 



the formation of the quarkonium state H from a cc pair in the state n. 



The leading order color-singlet processes are Z — > ipgg O] and Z — > ^cc [15|. Although 
both processes are of order a^, the latter is two orders of magnitude larger. The latter process 
has been interpreted as a fragmentation process in which the cc pair was first produced almost 
on-shell from the Z decay and then followed by the fragmentation of c or c into the if) [|16j . 
Therefore, the process Z — > ipcc is not suppressed by the quark propagator effect of order 
1/Mz as it does in the process Z — > ipgg. In the fragmentation approximation, the energy 
distribution of ip in the process Z — > ipcc is given by 

— {Z - if)(z) cc) ^2T{Z^ cc) x D c ^{z) , (3) 



D c _^) = sl n f K s 11 4- (16 - 32^ + 72z 2 - 32/ + 5z 4 ) , (4) 



where 



243m 3 (2 - z) f 

and z = 2E^/Mz with £^ denotes the energy of if). Both processes involve the same 
matrix element (Of( 3 Si)) whose value can be extracted from the leptonic width to be about 
0.73 GeV 3 . Numerically, the widths for Z — > ipgg and Z — > z/>cc are about 6 x 10~ 7 GeV 
and 7 x 10 -5 GeV, respectively. There is also a higher order color-singlet process Z — > ggg* 
followed by the gluon fragmentation g* — > -^gg [ITS . The branching ratio was estimated to 



be of order 10~ 6 only |5|, lSj . However, the recent results from LEP fl4j showed a branching 
ratio of order 10 -4 for prompt if) production, which is well above all the predictions from the 
color-singlet model. 

The leading order color-octet process is of order a s given by the process Z — > if>g, for 
which one of the Feynman diagrams is shown in Fig. 1(a). It is understood that soft hadrons 
are around to match the color. The decay rate of this channel is given by 

V17 h\ 647T 2 (0|(%)) ,. 167T M?S X )) 2 



2567T 2 (Q^( 3 Pq)) 
9 ^ M Z M% 



'■,••> p ,\ rq _ 3 £)2 6(1 + 2(1 + 3£ + 6£ 2 ^ 
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where £ = M|/M|, F = aa s /(x w (l — x w )), g v = l- — ~x w , g a = — |, and we have used a heavy 
quark spin symmetry relation (£^( 3 Pj)) « (2J + 1)(C^( 3 P )). This color-octet process has 
a very distinctive signature that the ip has an energy equal one-half of the Z mass and is 
recoiled by a hard gluon jet. Unfortunately, the short-distance factors are suppressed by at 
least one power of 1/Mz- Numerically, the width for Z — > ipg is of order 10 -7 GeV which 
renders this process useless at the Z resonance. 

The dominant color-octet process actually begins at the order a 2 in the process Z — > ipqq, 
for which one of the Feynman diagrams is shown in Fig. 1(b). Here q represents u, d, s, c, or 
b. The energy distribution of i[) for the process Z — > ipqq is calculated, in the limit m q = 0, 
to be 

dT , 



dz 



= ^nZ^qq)^^ 



1 



rrv c 

z z \z — Z L 



-2z L 



(6) 



z 2 — 4£) 1//2 . The physical range of z is 



where z = 2E i> /M z , £ = M|/M|, and z L = 
2y / £ < z < 1 + £. Similar expression was derived before for the process Z — > in QED 

P~5| , PU| . The above exact result can be simplified into the following form in the limit £ — ► 0: 



— (Z - ^(z)gg) « -^-^r(Z 



gg) 



(Of( 3 5x)> 



(z 



2z 



(7) 



Following the same procedures in Ref. [TO], one can show that the above limit corresponds 
to the sum of the fragmentation processes of Z — > qq followed by q — > ip or q — > ip and 
Z — > ggg followed by (7 — > ^. 

Using the value for the color-octet matrix element (Of ( 3 5i)) ~ 0.015 GeV 3 extracted 



from the fit to prompt ip data at the Tevatron [|TT| , |12l , m c = 1.5 GeV, »s 2m c , and 



a s (2m c ) = 0.253, we show the energy distribution of ip for this process in Fig. 2 (solid 
curve). The decay width of Z — > ipqq can also be obtained by integrating Eq. (|) over the 
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{ 5 (i_f)_2£log£ + 



2Li 2 [— ^— I - 2Li 2 



:i+o 2 } 



physical range of z, and is given by 

r(Z — > gg) 36 mj? 

-21og(l + log£ + 31og£ + log 2 £ 

~ 2.2 x 10" 4 . (9) 

Here Li 2 (x) = — Jq y log(l — £) is the Spence function. Summing over all the quark flavors 
(q = u,d,s,c,b), we obtain the decay width X) 9 r(Z — > ipqq) ~ 3.8 x 10~ 4 GeV and the 
branching ratio J2 q Br(Z — > ipqq) — 1.5 x 10~ 4 . Assuming the dominant prompt ip production 
process to be Z — > qqg* with g* — > ip + X, in which according to color-singlet model the 
off-shell gluon fragments into a ip plus two perturbative gluons, DELPHI obtained the limit 
Bt(Z — > qqg*; g* — >• ip + X) < 4.1 x 10~ 4 |T^]. Since the event topology of our color-octet 
process is similar to this one, this limit should also be valid for our color-octet process. 
Therefore, our result is consistent with this data. 

In Fig. 2 we compare the energy distributions of ip coming from the most important 
color-octet process Z — > ipqq using Eq. (j|) and the most dominant color-singlet process 
Z — > ipcc using Eq. For the color-octet process we have summed over all quark flavors, 
q = u,d, s, c, b, in the final state. The comparison in Fig. 2 shows a very spectacular difference 
between the color-octet and color-singlet processes. The spectrum for the color-octet process 
is very soft with a pronounced peak at the lower z end, while the spectrum for the color- 
singlet process is rather hard because of the nature of heavy quark fragmentation. Though 
the spectrum for the color-singlet process will be softened somewhat by including higher 



order QCD corrections [0, the spectacular difference remains unchanged. Another striking 
result is that the color-octet spectrum dominates over the color-singlet one for all values of 
z. This indicates that this color-octet contribution should not be neglected! 

The energy distributions for the ip' and \ c coming from the color-octet and color-singlet 
processes are similar to those illustrated above for the ip. The importance of the color- 
octet contributions can also be shown for these two cases as well. The rates for direct ip' 

7 



production are roughly scaled down by a factor of 5 from those of the if) shown in Fig. 2. 
Furthermore, the production of T by the color-singlet process Z — > Tbb and the color-octet 
process Z — > Tqq are similar to that of if), but with smaller branching ratios. Details will be 
presented elsewhere. 

In closing, we have considered the color-octet quarkonium production at the Z resonance. 
The dominant direct production of if) is shown to be the color-octet process Z — > ipqq, which 
is a few times larger than the most dominant color-singlet production process Z — > ipcc. 
Furthermore, the energy spectrum of the if) for the color-octet mechanism has been shown to 
be very different from that of the color-singlet process. Therefore, measurements of the decay 
width and the energy spectrum for the prompt if) production at LEP can further confirm the 
importance of the color-octet mechanism. 

This work was supported in part by the United States Department of Energy under Grant 
Numbers DE-FG02-84ER40173, DE-FG03-93ER40757, and DE-FG03-91ER40674. 

Note added: After this paper was submitted, similar results were reported by Cho |2TJ . 
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Figure Captions 

1. One of the contributing Feynman diagrams for the processes (a) Z — > ipg, (b) Z — > ipqq, 
and (c) Z — > ^cc. 

2. The energy spectrum dT/dz of the ^ from the color-octet process Z — > ^59 and the 
color-singlet process Z — > ipcc. 
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dr/dz (MeV) 




